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Abstract 

A  theoretical  analysis  is  presented  for  the  dynamic  behavior  of  a 
simply  supported  rigid,  perfectly  plastic  circular  plate  subjected  to  a 
rectangular  pressure  pulse.  It  is  shown  that  this  theory,  which  considers 
the  simultaneous  influence  of  membrane  forces  and  bending  moments ,  predicts 
final  deformations  which  are  considerably  smaller  than  those  given  by  the 
corresponding  bending  theory  of  Hopkins  and  Prager  even  when  maximum  de¬ 
flections  only  of  the  order  of  the  plate  thickness  are  permitted.  It  is 
believed  that  this  theoretical  analysis  should  assist  in  the  interpretation 
of  the  dynamic  biaxial  stress-strain  characteristics  of  materials  recorded 
on  diaphragms  fitted  in  impact  tubes  and  could  be  developed  further  in 
order  to  describe  the  behavior  of  plates  having  other  support  conditions 
and  different  dynamic  loading  characteristics. 
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plate  thickness 

impulse  per  unit  area  of  plate 
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<!uHpo) 


1/2 


o  H2/4 
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radial  and  circumferential  bending  moments  per 
unit  length 

o  H 
o 

radial  and  circumferential  membrane  forces  per 
unit  length 

transverse  shear  force  per  unit  length  of  plate 
outside  radius  of  plate 
principal  radii  of  curvature 


time  at  which  plate  reaches  permanent  position 

uniform  distributed  pressure  per  unit  area  of 
undeformed  plate 

dimensionless  bending  moments  M  /M  ,  M„/M 

6  r  o’  6  o 

dimensionless  membrane  forces  N  /N  ,  N„/N 

r  o’  0  o 

-k  sin$ 

6H 
_ o 

R2 

-k  cos<)) 

radial  coordinate  of  plate 
time 

displacement  in  direction  r  of  undeformed  plate 

transverse  deflection  perpendicular  to  undeformed 
plate 


Notation  (continued) 
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radial  and  circumferential  strains 

circumferential  coordinate  lying  in  plate 
radial  and  circumferential  curvatures 

k  /p 
o  *o 

mass  per  unit  area  of  plate 
radius  of  hinge  circle 

yield  stress  in  simple  tension 


t 


duration  of  pulse 


* 

(’) 

(  )' 


slope  of  the  mid-plane  of  a  plate  measured  in 
a  plane  which  passes  through  r  =  0  and  is 
perpendicular  to  the  plate  surface 


) 


1 .  Introduct ion 


The  behavior  of  rigid-plastic  circular  plates  under  the  influence  of 
static  loads  which  produce  infinitesimal  deflections  is  fairly  well  estab¬ 
lished  [1,  2,  3,  etc.].  When  finite  deflections  are  permitted,  however,  it 
is  observed  that  plates  can  support  external  loads  considerably  larger  than 
those  predicted  by  these  theories.  Onat  and  Haythornthwaite  [4]  indicated 
that  this  increased  load  carrying  capacity  is  due  mainly  to  the  important 
role  which  membrane  forces  play  in  the  finite  deformation  of  plates. 

It  is  clear  from  a  survey  of  the  pertinent  literature  that  most  atten¬ 
tion  has  been  directed  towards  the  dynamic  deformation  of  plates  in  which 
either  membrane  forces  [5,  6,  etc.]  or  bending  moments  [7,  8,  9,  etc.] 
alone  are  believed  to  be  important.  Moreover,  with  the  exception  of  some 
numerical  work  [10],  the  analysis  of  an  annular  plate  by  Florence  [11],  and 
some  recent  work  [12],  no  investigations  have  been  conducted  into  the  inter¬ 
action  effects  between  membrane  forces  and  bending  moments ,  although  such 
interaction  influences  considerably  the  static  loading  of  plates  [4]  and  the 
dynamic  loading  of  beams  [13].  Florence  [14]  applied  uniform  distributed 
impulses  to  some  simply  supported  circular  plates  and  observed  that  the  ap¬ 
propriate  rigid-plastic  theory  [8]  overestimated  considerably  the  recorded 
deflections  particularly  for  large  impulses.  Recently  it  has  been  demonstra¬ 
ted  [12]  that  a  significant  improvement  in  the  theoretical  predictions  of 
plates  loaded  impulsively  can  be  achieved  if  the  influence  of  membrane  forces 
and  bending  moments  is  retained  in  the  theory. 
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Symonds  [15]  indicated  that  the  permanent  deformation  of  rigid-plastic 
beams  subjected  to  central  force  pulses  having  rectangular  and  triangular 
shapes  differed  about  *15%  from  an  equivalent  half  sine  wave  pulse  with  the 
same  maximum  value  and  impulse.  Perzyna  [16]  developed  further  the  theory 
of  Hopkins  and  Prager  [7],  in  which  membrane  forces  are  disregarded,  and 
showed  that  for  a  given  impulse  the  character  of  the  pressure-time  function 
had  little  influence  on  the  final  shape  of  a  rigid-plastic  circular  plate. 
Hodge  [17]  and  Sankaranarayanan  [18],  on  the  other  hand,  found  that  the 
blast  characteristics  had  a  profound  effect  upon  the  final  deformation  of 
cylindrical  and  spherical  rigid-plastic  shells. 

In  practice,  the  blast  load  which  acts  on  a  plate  or  structure  often 
persists  for  a  short  period  of  time  rather  than  behaving  like  a  pure  im¬ 
pulse  as  assumed  in  Ref.  [12].  It  is  the  purpose  of  this  article,  there¬ 
fore,  to  study  the  behavior  of  a  rigid,  perfectly  plastic  circular  plate 
when  subjected  to  a  rectangular  pressure  pulse  such  as  the  one  shown  in 
Fig.  1.  The  results  of  this  analysis  will  be  compared  with  the  correspond¬ 
ing  values  from  Ref.  [7]  so  that  they  indicate  the  importance  of  membrane 
forces  and  with  those  of  Ref.  [12]  in  order  to  examine  and  assess  the  dif¬ 
ference  between  the  permanent  deflections  corresponding  to  a  pure  impulse 
and  an  equivalent  rectangular  pressure  pulse. 

2 .  Equilibrium  Equations 

It  may  be  shown  that  the  equilibrium  equations  for  the  finite  deflec¬ 
tions  of  a  circular  plate  subjected  to  axisymmetrical  dynamic  loads  can  be 
written  in  the  form  [12] 

(a0Nr)'  "  W  ara0Q/Rr  +  arV  +  ya0*  sin*  '  ya0a0,{i  cos*  = 


0 


(1) 
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(°eQ)’  +  arae^Nr/Rr  +  *VRe-l  +  araeq  +  naea  'w  cos<^  +  sin*  =  0  ^ 

<Vr>'  -  VMe  -  Ve°  =  0  <3> 

provided  the  rotary  inertia  effect  is  disregarded,  and 

a  =  1  +  e 
r  r 

aQ  =  r  +  u  =  r(l  +  eQ  ) 

1/Rp  =  *7(1  +  er) 

1/R0  =  sin<Ji/r 

The  positive  directions  of  the  various  quantities  are  indicated  in 
Fig.  2. 

If  we  limit  our  discussion  to  plates  having  small  strains  and  deflec¬ 
tions  which  are  not  too  large,  then  we  may  let  aQ  =  r,  =  1,  1/Rr  =  <f>'  , 
1/R0  =  sin<f>/r,  and  a0'  =  cos$  which,  using  cosiJ>  =  1,  and  sin<j>  =  -w'  , 
allow  equations  (l)-(3)  to  be  recast  as  follows 

rn  ’  +  n  -  n„  =  -rkw'/N  +  urww'/N  +  urii/N  (4) 

r  r  o  o  o  o 

and 


rm  "  +  2m'  -  m  '  -  4n  w'/H  =  rk/M  -  yrw/M  +  yruw'/M  (5) 

”  r  0  0  o  o  o 


where , 


n  =  N  /N 
r ,o  r , 0  o 


rn  ft/M 

r,0  r,0  o 


and  r$’Q,  rN^$'  and  <j>'w'  have  been  disregarded. 
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3.  Strains  and  Curvatures 


It  may  be  shown  for  small  strains  [19]  that 


e  =  (1  +  u' )u'  +  w'w' 
r 

h  =  “/r 

K  =  (1  +  u' )w"  +  u'w"  -  u"w'  -  u"w' 
r 

and 


=  w'/r 

0 


(6) 

(7) 

(8) 

(9) 


Griffith  and  Vanzant  [20]  observed  that  the  material  of  a  circular  plate 
tends  to  move  in  a  transverse  sense  at  high  rates  of  dynamic  loading.  This 
suggests  that 

u  =  0  (10) 


4.  Yield  Condition 

It  has  been  found  that  disregarding  elastic  effects  when  analyzing  canti¬ 
lever  beams  loaded  dynamically  is  a  powerful  simplification  and  a  valid  ap¬ 
proximation,  provided  the  external  energy  is  at  least  three  times  larger  than 
the  strain  energy  absorbed  by  the  beam  at  the  elastic  limit  [21].  Further, 
Frederick  [5]  and  Boyd  [6]  investigated  the  deformation  of  membranes  made 
from  work -hardening  material  and  found  that  a  simplified  perfectly  plastic 
analysis  provided  a  remarkably  accurate  model  of  the  true  behavior.  Conse¬ 
quently  the  plate  shown  in  Fig.  2  is  assumed  for  the  purposes  of  this  analysis 
to  be  made  from  a  rigid,  perfectly  plastic  material. 

The  yield  condition  proposed  by  Hodge  [22]  and  illustrated  in  Fig.  3  will 
be  used  in  this  article  since  it  simplifies  considerably  a  previous  analysis. 
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the  results  of  which  agree  reasonably  well  with  experimental  values  recorded 
on  plates  loaded  impulsively  [12],  This  approximate  yield  surface  is  an 
"upper"  bound  to  the  Tresca  yield  condition  for  a  uniform  shell  [23],  while 
a  similar  one  0.618  times  as  large  provides  a  "lower"  bound. 

5.  General  Equations 

Consider  a  rigid,  perfectly  plastic  circular  plate  which  is  simply  sup¬ 
ported  around  its  outer  edge  and  subjected  to  an  axisymmetrical  dynamic  load 
k(t)  ,  where  k(t)  is  a  function  of  time  and  is  transverse  to  the  mid-plane 
of  the  plate. 

Symmetry  demands  that  at  r  =  0 

m  =  m  =  -1  ,  and  n  -  n  =  1  ,  (11) 

r  o  t  u 

while  for  0  $  r  $  R  it  may  be  shown  that  equations  (6)-(10)  and  a  trans¬ 
verse  displacement  of  the  form 

w  =  W(t)(l  -  £)  (12) 

are  consistent  with  the  flow  rule  corresponding  to  the  yield  condition  illus¬ 
trated  in  Fig.  3  when 

m  =  -1  ,  n  =  1  ,  -1  $  m  $  0  and  0  $  n„  $  1  .  (13) 

0  r  r  8 

Substituting  equations  (10)  and  (13)  into  (4)  and  (5)  yields 

M  N  w' 

—  x—  (r2m  ')  =  k(t)  -  yw  +  — —  (14) 

2  r  r 

r4 

where  k(t)  is  an  external  load,  yw  is  an  inertia  term,  and  N  w'/r  is  in¬ 
troduced  when  finite  deflections  are  allowed.  It  may  be  shown  that  when  the 
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N^w'/r  term  is  disregarded  and  either  k(t)  or  pw  or  k(t)  -  pw  are 
retained,  then  equation  (14)  yields  the  same  results  as  quoted  in  refer¬ 
ences  [1,8,7],  respectively.  If  k(t)  and  N_w'/r  are  retained  and 
w  =  0  ,  then  equation  (14)  predicts  results  similar  to  Onat  and  Haythom- 
thwaite  [4]  for  deflections  at  r  =  0  greater  than  H/2  .  The  impulsive 
loading  case  in  reference  [12]  was  analyzed  using  equation  (14)  with 
k(t)  =  0  ,  while  this  article  is  concerned  with  dynamic  loading  for  which 
all  three  terms  must  be  included. 

Substituting  equation  (12)  into  (14)  and  integrating  gives 


where , 


„  T2„(t)  =  s  t 
a+2  V 


P 


(15) 


and  the  constants  of  integration  have  been  evaluated  using  the  conditions 

that  m  =  -1  at  r  =  0  ,  and  m  =0  at  r  =  R  . 
r  r 

It  may  be  shown,  using  equations  (12)  and  (14),  that 


„  _  k(t )  p_  d2W(t)  .1  _r. 

r  "  3M  M  ,2  o  ~  2R; 

o  o  dtz 


(16) 


which,  using  (15),  indicates  that  m  "  =  0  at  r  =  0  if 

r 


k(t )  _  2  + 


(17) 


Thus,  if  ^4^-  >  2  + 


‘  o 

4W(t  ) 
H 


,  then  the  yield  condition  given  by  equa¬ 


tion  (13)  will  be  violated  and  some  alternative  yield  condition  must  be  sought. 
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6 .  Rectangular  Pressure  Pulse  1  <  A  $  2 


The  rectangular  pressure  pulse  illustrated  in  Fig.  1  may  be  described 


viz 


and 


k(t)  =  k  for  0  $  t  5  T 
o 


k(t)  =  0  for  t  $  t  S  T 


(18) 

(19) 


where  r  is  the  duration  of  a  pressure  of  magnitude  kQ  and  T  is  the 
time  at  which  the  plate  finally  comes  to  rest. 

6.1  First  Stage  0  $  t  $  t 

The  general  solution  of  equation  (15)  for  this  case  is 

W(t)  =  A  cosrt  +  B  sinft  +  S/yz  +  2kQ/|iY2  (20) 

where  the  unknown  coefficients  A  and  B  may  be  determined  from  the  initial 
conditions 

w  =  w  =  0  at  t  =  0  . 

At  the  end  of  the  first  stage  equation  (20)  gives 


w  =  (X  -  1)(1  -  cosytKI  -  x) 


(21) 


and 


w  =  (A  -  DsinYT  (1  -  x) 


(22) 


where  x  =  r/R  ,  and  X  =  k  /p  . 

o  o 


6.2  Second  Stage  t  $  t  $  T 

The  general  solution  of  equation  (15)  for  this  stage  is 


W(t) 


C  coSYt  +  D  sinYt  +  6/y2 

TEOBimMS,  iimxur  ,  1 

BIOS.  013 

4BI8S1BM  FliOv'iui  sscuasj  ] 
SIEAP-H, 


(23) 
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where  the  unknown  constants  C  and  D  may  be  found  from  the  continuity 
requirements  which  demand  that  the  displacement  and  velocity  given  by  equa¬ 
tions  (12)  and  (23)  at  t  =  t  should  match  equations  (21)  and  (22). 

Thus , 


H 


w  =  —  [{(1  +  A(cosyT  -  1 ) }cosyt  +  Asinyrsinyt  -  1](1  -  x)  (24) 


but  w  =  0  at  t  =  T  . 
Hence , 


tanyT  =  mfcos^r  - 1) 


(25) 


and  at  t  =  T 


W  =  J  [/l  +  2A(1  -  cosyx ) ( A  -  1)  -  1](1  -  x) 


(26) 


which,  according  to  equations  (17)  and  (18),  is  valid  provided  1  $  A  $  2  . 

If  the  w'  term  is  disregarded  in  equation  (14),  then  it  may  be  shown 

that 

w  =  £  Ay2T2(A  -  1)(1  -  x)  ,  at  t  =  T  (27) 


which  is  the  same  as  the  result  obtained  by  Hopkins  and  Prager  [7]. 


7.  Rectangular  Pressure  Pulse  A  i  2 

It  has  been  shown  previously  that  the  yield  condition  given  by  equation 
(13)  is  not  suitable  for  rectangular  pressure  pulses  with  A  >  2  .  However, 
a  study  of  the  behavior  at  A  =  2  suggests  that  the  yield  condition  given 
by  equation  (11)  spreads  out  to  some  radius  r  =  pQ  ,  the  position  of  which 
is  related  to  A  .  The  behavior  of  a  thin  circular  plate  loaded  dynamically 
with  a  rectangular  pressure  pulse  having  A  z  2  can  be  considered  in  three 
stages ,  viz 
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1.  0$t$T.  k(t)  =  k  ,  and  a  stationary  hinge  is  formed 
at  r  =  p 

o 

2.  t  s  t  $  .  k(t)  =  0  ,  and  the  hinge  formed  in  part  1 

travels  inwards  from  r=p  at  t=T  to  r=0  at 


3.  $  t  $  .  k(t)  =  0  ,  and  the  hinge  remains  stationary 

at  r  =  0  until  the  plate  reaches  a  permanent  position  at 


7.1  First  Stage  0  $  t  $  t 


If  t  is  small  then  one  might  expect  that  the  w'  term  in  equation  (14) 
would  only  be  important  in  the  second  and  third  stages .  Therefore ,  disregard¬ 
ing  the  w'  term  in  equation  (14)  and  integrating  the  result  with  the  condi¬ 
tions  that  m^  =  0  at  r  =  R  and  m^  =  -1  at  r  =  pQ  ,  it  can  be  shown  that 


(28) 


if. 


kot2 

=  ir  >  f°r  0  s  r  5  Co 


(29) 


and 


w 


k  t2 

o  (R  -  r) 

2 M  (R  -  PQ ) 


for  p  $  r  $  R 
o 


(30) 


At  the  end  of  the  first  stage  t  =  x  ,  equations  (29)  and  (30)  give 


w  = 


k  t2 

o 

2u 


w  = 


k  x 
o 


,  for  0  $  r  $  p 


(31) 


and 


w 


V2  (R  -  r)  .  _  koT  (R  -  r) 

2y  (R  -  p  )  ’  w 


u  (R  -  pT  ’  for  Do  *  r  *  R  (32) 
o 
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7.2  Second  Stage  t  j:  t  £  T 


It  may  be  shown  that 


w  = 


k  t 
o 


for  0  $  r  s  p ( t ) 


(33) 


and 


*  -  k°T  (R  ~  r) 

W  u  (R  -  p(t))  ’ 


for  p(t)  S  r  $  R 


(34) 


are  consistent  with  the  flow  rule  for  the  yield  conditions 


m  =  m_=  -1  ,  n  =  n.  =  1  ,  when  0  $  r  S  p(t)  (35) 

r  o  r  0 

and 

m_  =  -1  ,  n  =1,  -Ism  *0,  0  $  n.  $  1  ,  when  p(t)  $  r  $  R 

0  r  r  0 

(36) 

Equations  (33)  and  (34)  are  continuous  across  r  =  p(t)  and  at  t  =  t 
match  the  values  given  by  equations  (31)  and  (32)  at  the  end  of  the  first 
stage . 

Substituting  equation  (34)  into  (14)  and  disregarding  the  w'  term, 
it  may  be  shown  that 

t  =  h  (1  - 1  -  (if  +  (s)3}  <37) 

which  suggests  that  an  appropriate  form  for  the  time  function  retaining  w' 
in  equation  (14)  is 

t  =  ’  ♦  T1  {1  - +(r)3  -  r  1  <38) 

Equation  (38)  gives  t  =  x  at  p  =  pq  and 

Ti  =  T  +  Ti  (1  -  f> 


when 


p  =  0 


(39) 
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where  is  an  unknown  constant,  the  value  of  which  will  be  determined 


later. 


Differentiating  equation  (38)  with  respect  to  time  gives 

1  T*  <R  +  3p)p 


R  -  p 


(40) 


R3 


In  order  to  analyze  the  behavior  of  the  plate  in  this  stage,  it  is  conven¬ 
ient  to  divide  it  into  three  separate  regions  0  $  r  £  p(t)  ,  p(t)  £  r  £  pQ  , 


and  p  £  r  S  R  . 
o 


7.2.1  0  £  r  £  p(t) 


The  deflection  w  using  equations  (31)  and  (33)  is 

k  T2  ^ 


W  = 


2y 


k  T 
o 


dt 


(41) 


from  which 


w’  =  0 


(42) 


Thus,  utilizing  equations  (33)  and  (42)  equation  (14)  reduces  to 


or 


_  (r2m, )  =  0 


m  =  -1 
r 


(43) 


which  is  consistent  with  equation  (35). 


7.2.2  p(t)  S  r  5  p. 


The  transverse  displacement  w  at  radius  r  is 


k 

w  =  ~ —  + 


2  v 


t(r) 


k  t 
o 


dt  t 


V  (R  -  r) 


t(r) 


y  (R  -  p ) 


dt 


(44) 
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where  from  (  38 ) , 


t(r)  =  t  + 


t;(i  -  (f)2  +  (?) 


3  2 


(45) 


Thus,  using  (40),  equation  (44)  yields 


k  tT’  .  5  9 

t  o  1/d  n  .  3p2  3r2. 

W  =  -  (Rp  -  Rr  +  -Jj - s— ) 

yR3  2  2 


(46) 


which,  when  substituted  with  the  derivative  of  (34)  into  (14),  yields 


m 


k  tR2 
o 


12MoT|(l+3y)(l-y)3 


(2x2  +  ^ - x3  - 6y2  +  4y3) 


1  (12xy  -  4x2  +  18y2x  -  3x3  +  _  12/  -  24y3 )  -  1 

(47) 

where,  x  =  r/R  ,  y  =  P/R  ,  and  the  constants  of  integration  have  been 

evaluated  from  the  requirement  that  m  and  m  '  are  continuous  across 

r  r 

r  =  p  . 

7.2.3  p  $  r  $  R 
o 


Now, 


and 


kox2(R  -  r)  f*  k  t ( R  -  r)dt 

W  =  2y(R  -  p  )  +  n(R  -  p) 

°  x 

-  k  t2  k  tT' 

2U(R  -  Pq)+  ~~  <2R>  +  -  2RP0  -  V> 


(48) 


(49) 


It  may  be  shown,  using  (34)  and  (49)  that  equation  (14)  gives 
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m  = 


k  tR2 
o 


12MoT|(l-y)3(l+3y) 


(2x2  -  x3  +  ^ - ^ - 6y2  +  4y3)  - 


X  X 


k  T2(x-y  )2  k  tT’  4y  3 

yHx(l-y  )  +  liT  (12xy  +  18xy2  '  12xyo  “  18xyo  “  — 


9y„ 


+  12y2  +  24y 3  +  -  12y2  -  24y3)  -  1 


o  x 


(50) 


where  the  constants  of  integration  have  been  evaluated  from  the  requirement 

that  m^  and  m^'  given  by  equations  (47)  and  (50)  are  continuous  at  r  =  pQ 

However ,  m  =0  at  x  =  1  .  Thus , 

’  r  ’ 


k  TR2(l-6y2+  8y3  -  3y4 )  k  x2(l-y  )  k  tT' 
O  J  O  O  OX 


12Mo(l-y)3(l+3y)T| 


yH 


6yH 


(12y  +  6y2  -  20y3  + 


9y4  -  12yQ  -  6y2  +  20y3  -  9y4)  -  1  =  0 


(51) 


which  putting  y  =  0  gives 


r  ti2  i  p  f  ip  =0 
11  2  1  3 


(52) 


where 


k  x 

C,  =  (-12y  -  6y2  +  20y3  -  9y4) 

1  6yH  Jo  jO  J  O  o 


(53) 


-k  T2 

c0  =  — (1  -  y  )  -  1 
2  yH  J  o 


(54) 


P  _  Ax 
3  ‘  2 


(55) 


Thus, 


Ti 


-C2  *  H  -  4C1C 

2C, 


(56) 


T1  =  T  +  Ti^  " 


and 


(57) 
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7.3  Third  Stage  ^  {  t  {  Tj 

The  transverse  displacement  given  by  equation  (12)  may  be  used  during 
this  stage  because  the  hinge  remains  stationary  at  r  =  0  .  However  since 
the  relations  for  w'  are  given  by  different  expressions  across  the  plate, 
then  it  is  necessary  to  divide  the  plate  into  two  sections  0  $  r  $  pQ  and 
PD  {  r  j  R  , 


7.3.1  0  $  r  $  p 


Using  equations  (12)  and  (46)  with  p  =  0  ,  it  may  be  shown  that 


w'  = 


-VT1 

uR3 


(*-  ♦  2£)  - 


(58) 


which  with  equation  (12)  allows  (14)  to  be  rewritten. 


!?<rV> 


JL 

M 


(r2 .  |i, .  (Ri2.  5|i) .  swt* 

dt*  R  Ur3h  2  RH 


or 


m  =  - 


2P, 


—  (2x2  -  x3)  - 
dt2 


k  tT' 
o  1 

6pH 


(4x2  +  3x3)  - 


2Wr 

RH 


-  1 


(59) 


(60) 


where  the  constants  of  integration  have  been  evaluated  from  the  condition 

that  m  =  -1  at  r  =  0  . 
r 


7.3.2  p  $  r  *  R 
o 


_  k  T2 

w'  =  -  ° _ 

2p(R  -  pQ) 


VTi 


(Rp 

R3  ° 


W(t) 

R 


(61) 


which,  when  substituted  into  equation  (14)  with  (12),  yields 


3_ 

3r 


(r2m_' )  =~g- 


d2w^t)  (r2 

dt2 


2k  T2r 
o 

p(R  -  pq)H 


4k  tT' 

O  1 

uR3H 


(Rp, 


4Wr 

RH 


(62) 
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from  which  it  may  be  shown  that 

«  3 


m  = 
r 


y  d2W  ,rz 
Mo  dt2  6 


12R 


3r  4Rr 


P  14  p  2  p  ^ 

O  yo  +  Mo  ^ 


3R 


k  t2  p  2 

yH(R  -  p  )  (r  +  “  "  2°o} 


2k  tT'  3p  2  p  2  ou  p  2  .  . 

- °_A  (Rp  +  — — )(r  +  ~  -  2p  )  -  H  (r  +  -f-  -  2p  )  -  (|  -  ~)T  +  G 

yR3H  02  r  o  RH  r  °  r  PD 

(63) 

where,  in  order  to  make  and  m  f  match  equation  (60)  and  its  derivative 

across  r  =  p  , 
o 


yp. 


F  =  - 


h2u  ,  P„  4k  tT'  _  3P  2WP  ‘ 

Ui  (I  -  _  °  1  P  3/R  __o) _ o_ 

2  3  4R  yR3H  °  8  m 


dt 


(64) 


and 


yp 


G  =  - 


12M 


o  dt 


h2u  p  ^k  tT'  D  p  2WP 

(2  -  -2.) _ 2_Jl  p  2(5.  +  _£) _ £ 

2  R  1  „o3l,  o  l6  +  8J  RH 


-  1 


yR3H 


(65) 


Now  if  mp  =  0  at  r  =  R  then  it  may  be  shown,  using  (63)-(65),  that 


d2W 

dt2 


+  y2W  =  6 


(66) 


the  general  solution  of  which  is 


W  =  M  cosyt  +  N  sinyt  +  6/y2 


(67) 


where 


(68) 


and 


6  -  - 


2P, 


2p  k  i' 
o  o 

y2H 


P  k  tT' 

(1  -  y  )  -  -  (12y  +  6y  2  -  20y  3  +  9y4) 

°  3y2H  °  °  (69) 


At  the  beginning  of  the  third  stage  t  =  T  ,  W  =  o  and  the  velocity 
given  by  equations  (12)  and  (67)  must  match  equation  (34)  at  the  end  of  the 
second  stage . 
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Hence , 


**6cosyT.  k  xsinyT. 

W(t)  =  -  I  -  +  — — — - J  cosyt  +  | 


MY 


'k  xcosyT  <SsinyT  . 

°  -  -  )  sinyt  +  6/y2 


PY 


(70) 


and 


j ScosyT  k  TsinyT  } 

—  *  - T— /slnTtt 


w  = 


fkoicosyT^  SsinyT^ 


cosyt  (71) 


But  w  =  0  at  t  =  T2  where , 


k  xcosyT,  6sinyT, 

o  1  1 


tanyT2  = 


6^  cosyT^  -  k^sinyT^ 
Y  — - - - 


(72) 


Thus , 


k  t2  k  xT' 
o  ,  o  1 

W  =  -2T+  —  (1 

,  M  k  2t2  1/2 

+  {—+(—  +  — - )  Hi 

y2  y4  p2y2 


2  x2  x3 
X  “  2  "2 


-  x)  ,  for  0  $  x  $  y 


(73) 


and 


k  t2(  1  -  x)  k  t(1  -  x)TI  3y2 

o  o  1  ,  J  o  \ 

- I -  (yo  +  — > 

o 


w  =  2y(l  -  y  )  + 


_2  k  2t2  1/2 

+  {—  +  (—  +  — - )  }(1  -  x)  ,  for  y  $  x  $  1 


Y2  Y4 


P2Y2 


(74) 


At  r  =  0  ,  equation  (73)  gives 


2t2  i/2 


w 


max 


k  t2  k  xT,  t  72  k  •‘-t 

4-  ♦  -2-i  ♦  4-  ♦  («i  +  — — — ) 

u  y  Y2  Y4  p2y2 


(75) 


If  the  w'  term  is  disregarded  in  (14),  then  it  may  be  shown  that 


>*  =  XHY2t2  ( 3X  -  2  )/16  ,  at  t  =  T,  (76) 

TTlclX  Z 


which  is  identical  to  the  result  obtained  by  Hopkins  and  Prager  [7]. 
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8.  Discussion 

The  yield  condition  [22]  indicated  in  Fig.  3  circumscribes  the  yield 
surface  of  a  uniform  shell  which  yields  according  to  the  Tresca  criterion 
[23].  A  solution  which  is  obtained  using  this  simplified  yield  condition 
is  termed  an  "upper"  bound  while  one  calculated  using  a  yield  surface  0.618 
times  as  large  and,  therefore,  lying  everywhere  inside  the  exact  yield  sur¬ 
face  is  referred  to  as  a  "lower"  bound.  It  is  assumed  that  these  "upper" 
and  "lower"  bounds  would  straddle  the  true  solution  based  on  an  exact  yield 
surface . 

It  is  clear  from  the  results  plotted  in  Fig.  4  that  membrane  forces  in¬ 
fluence  considerably  the  permanent  deformation  of  a  simply  supported  circu¬ 
lar  rigid-plastic  plate  loaded  dynamically  with  a  rectangular  pressure  pulse. 
Hopkins  and  Prager  [7]  predict  that  the  final  deformation  of  a  plate  in¬ 
creases  with  increase  of  (k  /p  )  for  a  given  impulse.  However,  inclusion 

o  o 

of  the  membrane  forces  gives  rise  to  a  trend  in  the  reverse  sense  which  in¬ 
dicates  that  membrane  forces  become  increasingly  important  with  increase  in 
A  . 

The  permanent  deformations  predicted  by  the  theory  presented  herein  for 
rectangular  pressure  pulses  of  various  magnitudes  are  compared  in  Fig.  5  with 
the  results  obtained  from  ref.  [12]  for  equivalent  impulses.  It  is  evident 
that  a  rectangular  pressure  pulse  with  A  =  12  predicts  results  similar  to 
those  presented  in  [12]  for  an  impulse,  while  for  smaller  values  of  A  , 
larger  deflections  are  predicted  at  a  given  magnitude  of  I'  .  It  is  worth 
noting  in  passing  that  the  curves  for  the  impulsive  case  plotted  in  Fig.  5 
bound  closely  some  experimental  results  recorded  by  Florence  [14]  on  Aluminium 
plates  with  R/H  =  16  . 
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9.  Conclusions 

A  theoretical  analysis  which  retains  the  influence  of  bending  moments 
and  membrane  forces  has  been  presented  for  a  simply  supported  rigid-plastic 
circular  plate  loaded  with  a  rectangular  pressure  pulse.  It  can  be  shown 
that  this  theoretical  analysis  predicts  final  deformations  which  are  con¬ 
siderably  smaller  than  those  given  by  the  bending  theory  of  Hopkins  and 
Prager  [7]  even  for  maximum  deflections  only  of  the  order  of  the  plate  thick¬ 
ness.  It  may  be  shown  that  a  rectangular  pressure  pulse  with  X  =  12  gives 
a  similar  maximum  permanent  deformation  to  an  equivalent  impulse,  while  for 
smaller  values  of  X  the  permanent  deflections  are  larger. 

It  is  thought  that  the  theoretical  analysis  presented  here  could  be 
developed  further  in  order  to  describe  the  behavior  of  plates  having  other 
support  conditions  and  different  characters  of  loading.  However,  it  is  be¬ 
lieved  that  some  estimate  of  the  influence  of  strain-rate  effects  should  be 
made  perhaps  in  a  manner  similar  to  those  of  Wierzbicki  [24]  or  Perrone  [25] 
who  disregarded  membrane  forces . 
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THE  FINAL  DEFORMATION  OF  A  SIMPLY 
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FIG.  5  COMPARISON  BETWEEN  THE  FINAL 
DEFORMATIONS  OF  A  RIGID  -  PLASTIC 
CIRCULAR  PLATE  SUBJECTED  TO  A 
RECTANGULAR  PRESSURE  PULSE  OR  AN 
IMPULSE. 


